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Abstract: Researching into changes in thermal growing season has been one of the most important 
scientific issues in studies of the impact of global climate change on terrestrial ecosystems. However, few 
studies investigated the differences under various definitions of thermal growing season and compared 
the trends of thermal growing season in different parts of China. Based on the daily mean air 
temperatures collected from 877 meteorological stations over northern China from 1961 to 2015, we 
investigated the variations of the thermal growing season parameters including the onset, ending and 
duration of the growing season using the methods of differential analysis, trend analysis, comparative 
analysis, and Kriging interpolation technique. Results indicate that the differences of the maximum values 
of those indices for the thermal growing season were significant, while they were insignificant for the 
mean values. For indices with the same length of the spells exceeding 5°C, frost criterion had a significant 
effect on the differences of the maximum values. The differences of the mean values between frost and 
non-frost indices were also slight, even smaller than those from the different lengths of the spells. 
Temporally, the starting date of the thermal growing season advanced by 10.0—11.0 days, while the ending 
dates delayed by 5.0—6.0 days during the period 1961—2015. Consequently, the duration of the thermal 
growing season was prolonged 15.0—16.0 days. Spatially, the advanced onset of the thermal growing 
season occurred in the southwestern, eastern, and northeastern parts of northern China, whereas the 
delayed ending of the thermal growing season appeared in the western part, and the length of the thermal 
growing season was prolonged significantly in the vast majority of northern China. The trend values of 
the thermal growing season were affected by altitude. The magnitude of the earlier onset of the thermal 
growing season decreased, and that of the later ending increased rapidly as the altitude increased, causing 
the magnitude of the prolonged growing season increased correspondingly. Comparing the applicability 
of selected indices and considering the impacts of frost on the definitions are important and necessary for 
determining the timing and length of the thermal growing season in northern China. 
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1 Introduction 


The global-mean surface temperature records show a linear warming trend of 0.85°C during the 
period 1880-2012, and this increase in surface temperature is projected to likely be 0.3°C-4.8°C 
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globally for 2081-2100 relative to 1986-2005 (IPCC, 2013). Global warming has shifted the 
timing of seasons with earlier spring, later autumn and earlier peak photosynthesis (Burrows et al., 
2011; Barichivich et al., 2012; Xu et al., 2016), and has extended the length of growing season 
(Piao et al., 2007; Wang et al., 2016). Changes in the vegetation growing season could not only 
affect the structure and function of ecosystems, but also lead to variations of agricultural 
production and carbon sequestration potential, which in turn may affect global and regional 
climate (Linderholm, 2006; Piao et al., 2007; Mozafari and Torki, 2015; Wang et al., 2016). 
Therefore, researching into the change of vegetation growing season has been one of the most 
important scientific issues in the studies of global climate change and terrestrial ecosystems 
(Mozafari and Torki, 2015; Yang et al., 2016). 

Thermal growing season is a period of time when plant growth can occur in theory, in other 
words, the constraints of temperature on plant growth are cancelled and growth will take place if 
other environmental needs are met (Carter, 1998). It is generally expressed as the number of days 
when air temperatures are higher than a predefined threshold (Shen et al., 2012; Guo et al., 2016). 
For example, Carter (1998) defined the start of growing season as daily mean temperature >5°C 
for >5 days and the end as the 10-day running mean when daily temperature <5°C. Dong et al. 
(2012) considered the thermal growing season as the period beginning on the last day of the first 
6-day spell with daily mean temperature >5°C after the last spring frost and the end on the first 
day of the first 10-day spell with mean temperature <5°C. Guo et al. (2016) took the start of 
growing season as the last day of the first 5-day spell with daily mean temperature >5°C, and the 
end as the last day of the first 5-day spell with mean temperature <5°C. 

Under the background of climate change, the changes in the mean temperature are anticipated 
to be mirrored in the changes of the growing season. However, the relationships between climate 
and thermal growing season are regionally different, and there is no universal method which is 
widely applicable to all regions to define the thermal growing season (Linderholm, 2006). 
Moreover, the indices of thermal growing season reflect the climate change in different ways, and 
even a small change in the indices definition may lead to distinctly heterogeneous results in the 
same regions. Although several scholars have analyzed the thermal growing season in the whole 
China or regions of China (Song et al., 2010; Dong et al., 2012; Yang et al., 2013; Guo et al., 
2016), few studies have conducted an in-depth comparison of the changes of the growing season 
in different parts of China with various definitions of indices (Song et al., 2010). Thus, the 
objectives of this study are to compare the different definitions of thermal growing season 
parameters (starting date, ending date, and the length of the thermal growing season) and to 
analyze their variations and trends in northern China, based on daily mean air temperature data 
from 877 stations during the period 1961-2015. 


2 Materials and methods 


2.1 Study area 


Northern China lies in the mid-latitudes of the Northern Hemisphere. This study was conducted in 
the regions north of the 33°N, including Xinjiang, Inner Mongolia, Heilongjiang, Jilin, Liaoning, 
Hebei, Beijing, Tianjin, Shandong, Shanxi, Ningxia, most of Henan, Shaanxi, Gansu and Qinghai, 
and some northern parts of Jiangsu, Anhui and Tibet (Fig. la). Strongly influenced by the East 
Asian monsoon, the study area is mostly characterized by arid and semi-arid climates, with the 
rainy season occurring in summer (June—August). Generally, the annual mean temperature ranges 
from —4°C to 14°C, and the annual precipitation varies from 25 to 1000 mm. There is a variety of 
complex topography, including mountains, basins, plateaus, hills and plains, and the regional 
landscape consists of various biomes, such as grasslands, farmlands, forests, deserts, shrubs and 
so on (Sun et al., 2015). 
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Fig. 1 Administrative division (a) and spatial distribution of selected meteorological stations (b) in northern China 


2.2 Data collection 


Daily mean air temperature data were used to calculate the growing season parameters (starting date, 
ending date, and the length of the growing season). These data came from the National 
Meteorological Information Center, China Meteorological Administration with the time span 
1951-2015. Totally, there were 1218 surface weather observation stations available in northern 
China. Given the higher missing rate of the original data and the less number of stations before 1961, 
we only used those observation data obtained during the period 1961—2015. Meanwhile, stations 
with average missing rates of over 1% were also excluded from our studies. Finally, 877 out of 1218 
meteorological stations were used, which had an average missing rate of less than 1% during the 
whole year and in spring and autumn over the past 55 years. The spatial distribution of selected 
stations was shown in Figure 1b. The missing data from a single day was filled in by the average of 
its adjacent days before and after, and the missing values of two or more than two consecutive days 
were filled with the simple linear regression method according to the observation data of 
neighboring stations which had the highest correlation coefficient with the former data and had no 
missing data during the same period. 


2.3 Methods 


We calculated a set of growing season indices, including the starting date of the growing season 
(SD, the ending date of the growing season (EJ), and the length of the growing season (LI), from 
daily mean temperature data based on different criteria as Table 1. In this study, four indices were 
used to determine the starting date and another five indices were used to determine the ending date 
of the growing season. Starting and ending indices were then combined to six different indices for 
the length of the growing season. Since the threshold of 5°C daily mean temperature has been 
widely accepted for determining the thermal growing season (Carter, 1998, Frich et al., 2002, 
Jones et al., 2002; Walther and Linderholm, 2006; Song et al., 2010; Dong et al., 2012; Guo et al., 
2016), here we also used it in all of the indices. The main differences between indices were 
presented in the lengths of the spells when daily mean temperature was over 5°C (starting date) or 
below 5°C (ending date). All indices that were used are listed in Table 1. 

In order to evaluate the strength of the relationship between indices, we compared the 
differences in both the mean and the maximum values for different combinations of indices. The 
difference between indices was firstly counted year-by-year at each station. Then, based on the 
55-year data sequence, we calculated the long-term average as the mean value at each station, and 
the maximum value at each station was simultaneously selected from the 55-year data sequence. 
Finally, the regional mean value was calculated using the station-averaged method and the 
regional maximum value was directly selected from the 877 stations and each combination had 
only one maximum value. 

The linear trends of time-series indices were also investigated. According to the latitude and 
longitude of stations, the Kriging interpolation technique was applied to generate the spatial 
patterns of the difference and trend for the thermal growing season, and the results were displayed 
by Golden Software Surfer. Given the sparse meteorological stations in some parts of northeastern 
China and most parts of western China, the interpolation results in these areas may not be able to 
correctly reflect the real variations to some extent. 
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Table 1 Definition of the indices of starting date, ending date and length of the thermal growing season in this 


study 
Parameter Index Definition 
=5d>5°C (SI 1) 5-day spell with Tmean remaining above 5°C 
Starting 
date of the >5d>5°C (SI 2) 6-day spell with Tmean remaining above 5°C 
proving =5d>5°C Fr (SI3) 5-day spell after the last Tmean spring frost with Tmean remaining 
season above 5°C 
(SD >5d>5°C Fr (SI 4) 6-day spell after the last Tmean Spring frost with Tmean remaining 
above 5°C 
=5d<5°C (EI 1) 5-day spell with Tmean remaining below 5°C 
Ending >5d<5°C (EI 2) 6-day spell with Tmean remaining below 5°C 
eie ai be =10d<5°C (EI 3) 10-day running mean of Tmean falling below 5°C 
growing 


season Fr or =5d<5°C (EI 4) First auinar Tmean-frost or 5-day spell with T mean 
(ED) remaining below 5C i 
Fr or >5d<5°C (EI 5) First aua imer Tmean-frost or 6-day spell with Tmean 
remaining below 5°C 


=5d>5°C | =5d<5°C (LI 1) Start (SI 1), End (EI 1) 


Length of >54>5°C | >5d<5°C (LI 2) Start (SI 2), End (EI 2) 

te =5d>5°C Fr | =10d<5°C (LI 3) Start (SI 3), End (EI 3) 

growing 

A >5d>5°C Fr | =10d<5°C (LI 4) Start (SI 4), End (EI 3) 
(LD) =5d>5°C Fr | Fr or =5d<5°C (LI 5) Start (SI 3), End (EI 4) 


>5d>5°C Fr | Fr or >5d<5°C (LI 6) 


Note: Tmean, mean air temperature; Fr, frost. 


Start (SI 4), End (EI 5) 


3 Results 


3.1 Differences for combinations of growing season indices 
3.1.1 Starting date of the growing season 


Time series of the mean starting date of the growing season are displayed in Figure 2. SI 1 
(=5d>5°C) had the earliest start, while SI 4 (>5d>5°C Fr) had the latest start of the growing 
season, with the difference of 4.4 days in the mean values during the period 1961—2015 (Table 2). 
Meanwhile, SI 1 had an earlier start of 3.2 days than SI 2 (>5d>5°C). Two indices including the 
frost criteria, i.e., SI 3 (]=5d>5°C Fr) and SI 4 had later start of the growing season than those not 
including the frost criteria (SI 1 and SI 2), with the differences of 1.8 and 1.2 days, respectively 
(Fig. 2). Therefore, the differences in the mean values between a frost index and a non-frost index 
are relatively little. Including the frost criterion or not in the definition can provide smaller 
differences than those from the different lengths of the spells exceeding 5°C. 

At 877 stations in northern China, the largest difference of the mean values for the starting date 
of the growing season was 12.4 days (Table 2) at Wushaoling station (3045 m a.s.1.) in Gansu. For 
the differences of the maximum values, the regional averages ranged from 17.5 to 28.5 days, and 
the maximum value were 64.0 days in pairs of SI 2—SI 1 and SI 4—SI 1, both occurred at Sheyang 
station in Jiangsu. For all the six combinations of the indices, the regional averaged difference 
was 2.2 days for the mean values, and 22.1 days for the maximum values over the past 55 years. 
Thus, the low difference in the mean values can only reflect the range of values to some extent, 
and the difference in the maximum values between indices was large. 


3.1.2 Ending date of the growing season 


For the five ending indices of the growing season, EI 3 (10d<5°C) had the earliest end, and EI 2 
(>5d<5°C) had the latest end of the growing season, with the difference of 9.4 days for the mean 


JOURNAL OF ARID LAND 


= 
tes) 


~ | (b) = >5d>5°C (SI 2) 
96 j i we ech cl 


89 


Starting date (Julian day) 


—o L 4 1 4 L 4 L Jt 


68 n 4 4 n n L 1 4 68 
1961 1967 1973 1979 1985 1991 1997 2003 2009 2015 1961 1967 1973 1979 1985 1991 1997 2003 2009 2015 


Fig. 2 Effects of consecutive days and frost (Fr) criterion on the starting date of the growing season during the 
period 1961-2015 in northern China 


Table 2 Differences of the mean and maximum values for combinations of starting indices at stations in 
northern China 


Differences of the mean values (days) Differences of the maximum values (days) 
Combinations 
Mean Maximum Mean Maximum 

SI 2-SI 1 32) 1E 26.3 64 
SI3-SI 1 1.8 10.0 21.1 58 
SI4-SI 1 4.4 12.4 28.5 64 
SI 2-SI3 1.4 T2 21.6 (-17.5) 53 (-57) 
SI 4-SI 2 1.2 8.4 17.5 57 
SI 4-S1 3 2.6 We 21.6 53 


Note: For the mean differences, the positive and the negative values are counted together; and for the maximum differences, they are 
counted separately. The negative values in parentheses represent the former index is less than the latter index at some stations or in some 
years. The description of the difference is the same as in Tables 3 and 4. 


values during the period 1961-2015 (Table 3). Meanwhile, EI 2 had a later end of 1.8 days than 
EI 1 (=5d<5°C), and EI 1 had a later end of 7.6 days than EI 3 averagely. The difference between 
EI 4 (Fr or =5d<5°C) and EI 5 (Fr or >5d<5°C) was relatively low (1.1 days). Two indices 
including the frost criteria, i.e., EI 5 and EI 4 had earlier end of the growing season than EI 2 and 
EI 1, with the differences of 1.3 and 2.0 days for the mean values, respectively (Fig. 3). Therefore, 
indices including the frost criterion or not can also cause smaller differences than those from the 
different lengths of the spells in the definition of the ending date of the growing season. 

In northern China, the largest difference of the mean values for the ending date of the growing 
season was 14.1 days, at Hua Shan Mountain station in Shaanxi for the pair of EI 2 vs. EI 3 
during the period 1961—2015 (Table 3). For the differences of the maximum values, the regional 
averages ranged from 7.9 to 29.9 days, and the maximum value was 54.0 days occurred at 
Suiping station in Henan. For all the ten combinations of the indices, the regional averaged 
difference was 4.0 days for the mean values and 23.0 days for the maximum values, respectively, 
while for the six combinations excluding EI 3, the regional averaged differences for the mean and 
maximum values were 1.6 and 20.7 days. Thus, the difference among indices was in general small 
for the mean values, especially when EI 3 index was excluded, but that for the maximum values 
was still large. 

3.1.3 Length of the growing season 

For the six indices of the length of the growing season, LI 4 (>5d>5°C Fr | 10d<5°C) had the 
shortest length, and LI 1 (=5d>5°C | =5d<5°C) had the longest length of the growing season, with 
the difference of 12.0 days for the mean values (Table 4). Meanwhile, LI 1 had a longer growing 
season of 9.3 days than LI 3 (=5d>5°C Fr | 10d<5°C), and LI 2 (>5d>5°C | >5d<5°C) had a 
longer growing season of 10.5 days than LI 4 over the past 55 years. The difference for the pair of 
LI 1 vs. LI 2 was the least (1.4 days). Two indices, including the frost criteria, LI 5 (=5d>5°C Fr | 
Fr or =5d<5°C) and LI 6 (>5d>5°C Fr | Fr or >S5d<5°C), had shorter length of the growing season 
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than those not including the frost criteria (LI 1 and LI 2), with the differences of 3.1 and 3.2 days, 
respectively (Fig. 3). Once again it showed that including the frost criterion or not could cause 
smaller differences of the mean values than those from the different lengths of the spells in the 
definition of the length of the growing season. 


Table 3 Differences of the mean and maximum values for combinations of ending indices at stations in northern 
China 


eer Differences of the mean values (days) Differences of the maximum values (days) 
Mean Maximum Mean Maximum 

EI 2-EI 1 1.8 4.8 21.6 43 
EI 1-E1 3 7.6 11.6 26.6 54 
EI 1-EI 4 1.3 8.4 18.7 47 
EI 1-EI5 0.2 12.8 18.7 (—18.7) 47 (-43) 
EI 2-EI 3 9.4 14.1 29.9 54 
EI 2-EI 4 3.1 10.4 25.1 53 
EI 2-EI 5 2.0 9.8 PAIS) 53 
EI 4-EI 3 6.2 9.6 23.5 (10.6) 54 (-35) 
EI 5—EI 3 T3 12.7 26.2 (-7.9) 54 (-35) 
EI 5—EI 4 1.1 3.1 18.7 43 
321 321 
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Fig. 3 Effects of consecutive days and frost (Fr) criterion on the ending date of the growing season during the 
period 1961-2015 in northern China 


At 877 stations in northern China, the largest difference of the growing season's length for the 
mean values was 22.5 days at Wushaoling station in Gansu during 1961—2015, which occurred in 
the pair of LI 1 vs. LI 4 (Table 4). For the differences of the maximum values, the regional 
averages were ranged from 3.3 to 39.2 days, and the maximum value was 93.0 days occurred at 
Xifeng station (1421 m a.s.l.) in Gansu. For the whole study area, the regional averaged 
difference was 5.7 days for the mean values and 28.6 days for the maximum values, if all fifteen 
combinations were considered. Thus, although the average difference of the mean values was 
small, the difference in the maximum values was also large. 
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Fig. 4 Effects of consecutive days and frost (Fr) criterion on the length of the growing season during the period 
1961-2015 in northern China 
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Table 4 Differences of the mean and maximum values for combinations of length indices at stations in northern 
China 


Differences of the mean values (days) Differences of the maximum values (days) 


Combinations 
Mean Maximum Mean Maximum 
LI 1-LI 2 1.4 55 25.5 (19.8) 64 (—42) 
LI 1-LI3 95 175 33.4 78 
LI 1-LI 4 12.0 25 39.2 (-3.3) 85 (32) 
LI 1-LI5 3.1 14.3 26.7 73 
LI 1-LI 6 4.6 21.6 31.9 (-16.2) 73 (-42) 
LI 2-LI3 VS 21.0 33.7 (-15.1) 93 (-45) 
LI 2-LI 4 10.5 17.0 34.7 93 
LI 2-LI 5 1.7 U7) 28.5 (—20.5) 69 (—53) 
LI 2-L1 6 32 12.9 26.8 66 
LI 3-LI 4 2.6 2 21.6 53 
LI 5-LI 3 6.2 9.6 23.5 (—10.6) 54 (-35) 
LI 6-LI 3 4.7 9.0 24.5 (-18.2) 53 (-44) 
LI 5-LI 4 8.9 14.6 31.4 8.8) 71 (35) 
LI 6-LI 4 Wd 7 26.2 (-7.9) 54 (-35) 
LI 5-LI 6 1.5 4.5 21.2 (-17.0) 53 (-42) 


3.2 Spatial and temporal trends for growing season indices 


3.2.1 Starting date of the growing season 


Except for several scattered stations in western China, earlier onset of the growing season 

occurred at all stations, regardless of the indices selected. The difference of trends among four 

indices was relatively small, so the trend values were little affected by the method chosen (Fig. 5). 

For the four indices of the growing season, the starting date advanced significantly at rates of 

1.0-5.5 days per decade in the southwestern, eastern and northeastern parts of the study area, 

including western Qinghai, most of Inner Mongolia and Shaanxi, north-central Heilongjiang, 

western Jilin, southern Liaoning, southern Hebei, Shandong, northeastern Henan, northern Anhui, 

northern Jiangsu, western Shanxi, Ningxia, southeastern and north-central Gansu, and some areas 

in Xinjiang and Tibet, of which areas with earlier onset of the growing season by 1.0-3.7 days per 

decade accounted for the largest proportion. 

Variation of starting date 
(days per decade) 
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Fig. 5 Linear trends (days per decade) in the starting indices of the growing season during the period 1961-2015 
in northern China. Areas of significant changes (P<0.05) were colored while those with insignificant changes 
(P>0.05) were not colored and shown as white spaces. The description is the same as in Figures 6 and 7. 


~ NINAX IV PROT 


CUI Linli et al.: Variations of the thermal growing season during the period 1961-2015 in northern... 


For all 877 stations, the averaged changes in four starting indices of the growing season ranged 
from —11.3 to -10.1 days during 1961-2015 (Table 5). Statistics showed that there were 608, 602, 
586 and 576 stations with significant advancing rates of 1.0—5.5 days per decade for SI 1, SI 2, SI 
3 and SI 4, accounting for 69.3%, 68.6%, 66.8% and 65.7% of the total stations across northern 
China, respectively, of which there were 54.4%, 56.0%, 56.6% and 56.3% of stations with 
significant advancing rates of 1.0-3.0 days per decade, respectively. Areas with significant 
advance of the growing season reduced gradually from SI 1 to SI 4 (Fig. 5). None of station had 
significant later onset of the growing season for all four indices. Indices including the frost 
criteria brought about changes of 1.0 day less than those disregarding frost criteria on average. 


Table 5 Changes in the regional averaged growing season indices during the period 1961-2015 in northern 
China 


Average Change with altitude (days) 
Parameter Index change 
(days) 0-1000m 1000-2500 m 2500-4000 m 
=5d>5°C (SI 1) -11.3 -11.5 -10.9 -8.6 
gamadi CECEO -11.0 i2 -10.6 -8.4 
RONAN SERIO =5d>5°C Fr (SI 3) —10.2 -10.8 -8.9 = 
>5d>5°C Fr (SI 4) -10.1 —10.6 -8.8 -7.9 
=5d<5°C (EI 1) 5.1 4.5 6.3 10.2 
>5d<5°C (EI 2) 5.3 4.7 6.6 9.6 
SUAS UNS KE) 5.6 48 71 10.1 
growing season 
Fr or =5d<5°C (EI 4) ST 4.9 7.0 12.6 
Fr or >5d<5°C (EI 5) 59 5.1 7.4 12.1 
=5d>5°C | =5d<5°C (LI 1) 16.4 16.0 172 18.8 
>5d>5°C | >5d<5°C (LI 2) 16.3 15.9 17.2 18.0 
Lanke =5d>5°C Fr | =10d<5°C (LI 3) 15.8 15.7 15.9 17.8 
HONS SEESON >5d>5°C Fr | =10d<5°C (LI 4) Sa 15.5 15.9 18.0 
=5d>5°C Fr | Fr or =5d<5°C (LI 5) 16.0 15.8 15.9 20.2 
>5d>5°C Fr | Fr or >5d<5°C (LI 6) 16.0 15.8 16.2 20.0 


Over the past 55 years, the maximum difference for the change of the starting date of the 
growing season was 19.0 days at Shangnan station in Shaanxi among four indices used, and 
Fuping station in Hebei seemed to be the most insensitive to the indices chosen. The earlier 
starting date of the growing season was less in higher altitude areas than in lower areas for the 
four indices (Table 5). In areas with altitudes ranged from 0 to 1000 m a.s.1., the starting date of 
the growing season had averagely advanced about 10.6—-11.5 days, but in areas with altitudes 
ranged from 2500 to 4000 m a.s.l., it advanced about 7.7-8.6 days during 1961—2015. At Tongde 
station (3148 m a.s.l.) in Qinghai and Huining station (1739 m a.s.I.) in Gansu, the earlier onset of 
the growing season was the most remarkable (more than 27.0 days). 


3.2.2 Ending date of the growing season 


Later end of the growing season was observed in the vast majority of northern China except for 
some small areas in southern and western Liaoning, eastern and southern Hebei, eastern Shanxi, 
western and southern Henan, the north end of northeastern China and two stations (Aral station 
and Hoxud station) in Xinjiang. The difference of trends among five indices was also small, and 
the trends were little affected by the indices chosen (Fig. 6). Being different from the earlier onset 
of the growing season, which was significantly concentrated in the eastern and northeastern parts, 
the later end of the growing season was significantly concentrated in the western part of the study 
area. In northern Tibet, most of Xinjiang and Qinghai, Jilin, central and southwestern 
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Heilongjiang, and some areas in Inner Mongolia, Gansu, Shaanxi, Shanxi, Shandong and northern 
Jiangsu, the ending date was delayed significantly at speeds of 0.5-3.8 days per decade during 
1961-2015, of which areas with later end of the growing season by 1.0-3.0 days per decade 
accounted for the largest proportion. 


Variation of ending date 
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Fig. 6 Linear trends (days per decade) in the ending indices of the growing season during the period 1961-2015 
in northern China 


Generally, the change of the ending date of the growing season was small. For the five indices, 
the averaged changes ranged from 5.1 to 5.9 days over northern China during 1961-2015 (Table 
5). There were 297, 308, 430, 319 and 331 stations with significant delaying rates of 0.5-3.5 days 
per decade for EI 1, EI 2, EI 3, EI 4 and EI 5, accounting for 33.9%, 35.1%, 49.0%, 36.4% and 
37.7% of the total stations, respectively, of which there were 31.9%, 33.5%, 46.8%, 35.0% and 
36.3% of stations with significant delaying speeds of 1.0-2.5 days per decade, respectively. For 
the five indices, EI 3 had the largest areas with significant delay of the growing season, and EI 1 
had the smallest areas (Fig. 6). None of station had significant earlier end of the growing season 
for all five indices. Indices including the frost criteria produced changes of about 0.6 days more 
than those disregarding frost produced averagely. 

During the period 1961—2015, the maximum difference for the change of the ending date of the 
growing season was 25.9 days at the Wushaoling station among five indices, and Nomhon station 
(2790 m a.s.l.) in Gansu was the most insensitive to the indices chosen. The ending date of the 
growing season was delayed more in higher altitude areas than in lower altitude areas for the five 
indices (Table 5). In areas with altitudes ranged from 0 to 1000 m a.s.l., the ending date of the 
growing season had averagely delayed about 4.5-5.1 days, but in areas with altitudes ranged from 
2500 to 4000 m a.s.l., it averagely delayed about 9.6—12.6 days over the past 55 years. At Tongde 
station and Wushaoling station, the later end of the growing season was the most obvious (more 
than 20.0 days). 

3.2.3 Length of the growing season 


Except for several stations in Xinjiang, Heilongjiang and Hebei, extended length of the growing 
season occurred in the whole study area, regardless of the indices selected, and the difference of 
trends among six indices was also small (Fig. 7). In general, the length of the growing season was 
significantly extended at rates of 0—9.1 days per decade in the vast majority of northern China 
during the period 1961-2015, of which areas with prolonged length of the growing season by 
1.5—6.0 days per decade accounted for the largest proportion. 

The averaged changes for the six indices of the length of the growing season ranged from 15.7 
to 16.4 days during the period 1961—2015 (Table 5), and at most of the stations, they were 
extended from 8.0 to 25.0 days. There were 692, 697, 756, 749, 713 and 712 stations with 
significant prolonging rates of 0-9.1 days per decade for LI 1, LI 2, LI 3, LI 4, LI 5 and LI 6, 
accounting for 78.9%, 79.5%, 86.2%, 85.4%, 81.3% and 81.2% of the total stations across 
northern China, respectively, of which there were 70.2%, 69.4%, 80.2%, 78.8%, 73.3% and 


201803.00004v1 


chinaXiv 


ChinaXiv@ (ESAT! 


CUI Linli et al.: Variations of the thermal growing season during the period 1961-2015 in northern... 


71.2% of stations with significant prolonging rates of 1.5-4.5 days per decade, respectively. For 
the six indices, LI 3 had the largest areas with significant extension of the growing season, and LI 
1 had the smallest areas. Indices including the frost criteria caused the changes of 0.3—0.4 days 
shorter of the length of growing season than those ignoring frost criteria produced averagely. 

Over the past 55 years, the maximum difference for the change of the length of growing season 
was 28.1 days at the Wushaoling station among six indices used, and Dingxin station (1177 m 
a.s.l.) in Gansu was the most insensitive to the indices chosen. The length of the growing season 
was prolonged more in higher altitude areas than in lower areas for the six indices (Table 5). In 
areas with altitudes ranged from O to 1000 m a.s.l., the length of the growing season had 
averagely extended about 15.5—16.0 days, but in areas with altitudes ranged from 2500 to 4000 m 
a.s.l., it averagely extended about 17.8—20.2 days during the period 1961-2015. Again at Tongde 
station and Wushaoling station, the extension of the growing season was the most prominent 
(more than 40.0 days), and at Aral station (1012 m a.s.l.) in Xinjiang there were small decreases 
of the length of growing season for all the five indices (not significant). 


Variation of length 


80°E 90°E 100°E 110°E 120°E 130°E 80°E 90°E 100°E 110°E 120°E 130°E (days per decade) 
L 1 4 4 l 4 1 1 L 1 1 1 

53°N 53°N 53°N 53°N i 
(a) LI 1 TASSA NT ()LI2 S 3 15 
48°N | 48°N 48°N 6.0 
43°N 4 43°N 43°N ; 4.5 
38°N 38°N 38°N : ` 38°N [30 
2 1.5 
33°N dat 3 aa 
80°E 90°E 100°E 110°E 120°E 130°E 7 


80°E 90°E 100°E 110°E 120°E 130°E 
80°E 90°E 100°E 110°E 120°E 130°E 80°E 90°E 100°E 110°E 120°E 130°E 
53°N 53°N 


48°N 48°N 


43°N 43°N 
38°N 38°N 


33°N 33°N a 
80°E 90°E 100°E 110°E 120°E 130°E 80°E 90°E 100°E 110°E 120°E 130°E 


53°N 53°N 
48°N 48°N 


43°N 43°N 
38°N 38°N 
33°N OS AN. 1000 kn} 33°N 330N S 

80°E 90°E 100°E 110°E 120°E 130°E 80°E 90°E 100°E 110°E 120°E 130°E 


Fig. 7 Linear trends (days per decade) in the length of the growing season during the period 1961-2015 in 
northern China 


4 Discussion 


The length of growing season is important for local agriculture, carbon sequestration potentials, 
ecological environment and human activities (Brown et al., 2010). Given that the difference of 
mean values between a frost index and a non-frost index was relatively small, it had less influence 
on the statistical analysis of the growing season in northern China (Figs. 2—4) whether including 
the frost criterion or not. However, the difference of maximum values between a frost index and a 
non-frost index was large, so in this context the frost criterion had obvious effects on the 
determination of annual starting date, ending date and the length of the growing season at stations. 
This is consistent with existing results. For example, Walther and Linderholm (2006) noted that 
including or excluding the frost criteria had a significant impact on the starting date of the 
growing season in the south western parts of the Greater Baltic Area, though disregarding the 
frost criteria during autumn did not have the same importance for the ending date of the growing 
season. The results from Guo et al. (2013) also showed that frost criterion had a significant 
influence on the growing season in most central parts of Inner Mongolia, particularly on the onset 
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of the growing season. 

During the period 1961—2015, the starting date of the growing season had shifted 10.0—11.0 
days earlier while the ending date had moved 5.0—-6.0 days later, extending the length of the 
growing season by 15.0—16.0 days in northern China depending on the indices chosen (Table 5). 
These results are in agreement with others. Song et al. (2010) showed that the thermal growing 
season had extended at a rate of 2.3 days per decade in northern China during the period 
1951-2007. In the Inner Mongolia, China, Guo et al. (2014) also found that the length of the 
growing season had extended by 22, 14, 17, and 15 days in the desert, grassland, agriculture, and 
forest regions during the period 1961—2010, respectively. In temperate China, Shen et al. (2012) 
found that the starting date of the thermal growing season had shifted 8.4 days earlier though the 
ending date had moved 5.7 days later from 1960 to 2009 averagely. The results from Guo et al. 
(2016) also showed that from 1961 to 2013, there were 9.0 and 11.1 days earlier of the onset of 
tree growing season, and 4.8 and 5.8 days later of the ending date of tree growing season, 
extending the length of the tree growing season by 13.3 days in the mountainous regions of Hebei 
and Shanxi and 16.4 days in the Loess Plateau. 

In northern China, the impact of different indices from the same category on the long-term 
trend of the growing season was relatively small, but there were still some differences in different 
areas and the differences between indices were not homogeneous over the whole area (Figs. 5—7). 
Walther and Linderholm (2006) compared a lot of definitions of the growing season in the Greater 
Baltic Area, and also found some differences in trends of the length of growing season depending 
on the definition used. In addition, the 877 meteorological stations were not evenly distributed, 
and most of them located in the southeastern part of northern China (Fig. 1). In the southwestern 
part, including southern Xinjiang, northern Tibet, and southwestern Qinghai, the distribution of 
stations was sparse. Given that the most widespread land types are deserts, mountains and bare 
lands in these areas, we calculated the mean and maximum values in this paper using data from 
the meteorological stations, the relatively small number of stations available in these areas had 
less impacts on our results, though the interpolation of the results in these areas might have some 
uncertainties. 

Thermal growing season can only estimate the real growing season to a certain extent, and it 
should be seen as a model approach to determine the growing season using daily mean 
temperature (Walther and Linderholm, 2006). Phenological data obtained from ground 
observation can provide real and exact results of growing season parameters because it is the 
integration of all related environmental factors that promote the growth of plants in a particular 
area. Zheng et al. (2002) analyzed the changes in the phenology of plants in China based on 
phenological observation data from 26 stations and found spring phenophases had been advanced 
in northern China, northeastern China and the lower Yangtze region, but delayed in the middle 
reaches of Yangtze River and the eastern part of southwestern China since 1980s. Hao et al. (2017) 
analyzed the phenology data of typical ligneous plants in Kashgar Prefecture and results showed 
that the phenology of ligneous plants had shift to an earlier date in spring, on average, it had been 
advanced 8 days gradually during the period 1985—2012. Based on the phenological data of Larix 
gmelinii during 1987-2012, Yang et al. (2016) found that that growing season of L. gmelinii was 
prolonged in the Greater Khingan Mountains, especially the first leaf unfolding and the end of 
leaf fall were delayed by 17.3 and 12.1 days per decade, respectively. However, the availability of 
phenological data is more limited and only in shorter time scales. 

With the development of satellites, remote sensing data can now provide continuous coverage 
and fine spatial resolution of regional or global vegetation and thus widely used in phenology 
monitoring (White et al., 2005; Wang et al., 2016). For example, Hggda et al. (2013) analyzed the 
trends of growing season in Fennoscandia using the Global Inventory Modeling and Mapping 
Studies NDVI3g data and results showed an earlier start of the growing season of 9.8—13.8 days 
during the period 1982-2011. Nagai et al. (2015) examined the variations in the timing of the 
growing season in Japan using daily satellite-observed green-red vegetation index, and found the 
sensitivity of the ending date of the growing season to air temperature was much less than that of 
starting date in both deciduous broadleaf and needle leaf forests during the period 2003-2012. 
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Chi et al. (2016) showed the heterogeneous trends of the onset of the growing season for different 
vegetation types in Xilin Gol League, China by using NOAA/AVHRR NDVI data of 1-km 
resolution during the period 1989-2009. Wang et al. (2016) examined the onset, length, and 
ending date of the growing season in the Loess Plateau based on NDVI data, and results showed 
that 54.84% of the vegetation had an earlier starting date of the growing season, whereas 67.64% 
had a later end during the period 2000-2010. Therefore, it is the focus of future research to 
carefully select the indices and compare the applicability of different indices with a variety of 
temperature thresholds, combining satellite remote sensing monitoring and plant phenological 
observation. 


5 Conclusions 


In northern China, the differences of the maximum values among different parameters of the 
thermal growing season (i.e., starting date, ending date and the length of the growing season) 
were large, though the differences of the mean values between these indices were small. 
Therefore, it is necessary to carefully select the indices of the growing season and to compare the 
applicability of different indices in a variety of geographic and climatic regions. Indices including 
the frost criterion or not could cause little differences for the mean values, even smaller than those 
from the different lengths of the spells, but for the maximum values, the differences between frost 
and non-frost indices were large. Thus the consideration of frost criterion is important for 
determining the annual timing and the length of growing season in northern China. 

Over the past 55 years, the regional average trends of the thermal growing season showed an 
earlier onset of 10-11 days, a delayed end of 5.0-6.0 days and consequently a prolonged length of 
15.0—16.0 days in northern China. Spatially, the starting date of the growing season advanced 
significantly at rates of 1.0-3.7 days per decade in the southwestern, eastern and northeastern 
parts of northern China, whereas the ending date delayed significant at rates of 1.0—3.0 days per 
decade in the western part, and the length of the growing season was prolonged significantly at 
rates of 1.5—6.0 days per decade in the vast majority of northern China. As the altitude increased, 
the advancing rate of the growing season decreased, and the delaying rate of the ending date of 
the growing season increased rapidly, correspondingly, the prolonging rate of the length of the 
growing season increased. 
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